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Abstract —In networked systems, state estimation is hampered by 
commnnication limits. Past approaches, which consider schednling sen¬ 
sors through deterministic event-triggers, reduce communication and 
maintain estimation quality. However, these approaches destroy the 
Gaussian property of the state, making it computationally intractable 
to obtain an exact minimum mean squared error estimate. We propose 
a stochastic event-triggered sensor schedule for state estimation which 
preserves the Gaussianity of the system, extending previous results from 
the single-sensor to the multi-sensor case. 

I. Introduction 

Networked Control Systems (NCSs), spatially distributed systems 
where sensors, actuators, and controllers exchange information over 
a shared, bandlimited communication network, have become a topic 
of significant interest in both academia and industry. As noted by 
fUl . the use of NCSs in practice provides for flexible architecture 
and reduces costs in installation and maintenance. Thus, NCSs have 
been used in several applications including public transportation, 
health care, and mobile sensor networks. Nonetheless, remote state 
estimation remains a significant challenge in NCSs Q]. Traditionally, 
state estimates are computed at an estimation center using information 
from sensors which sample and send measurements periodically. 
While it is reasonable to assume that remote state estimation centers 
are well equipped, in most cases, sensors have a limited power supply 
and are difficult to replace. Moreover, bandwidth constraints in a 
communication network may restrict the number of sensors which 
can communicate at any given time 0, 0, 0. One way to address 
these issues is to simply reduce the communication rate. This solution 
however degrades estimation quality. In this paper, we propose a 
sensor scheduling scheme which allows us to achieve a desired 
tradeoff between communication rate and estimation performance. 
Specifically, we design a stochastic multi-sensor event-based schedule 
for the remote state estimation problem which extends the single 
sensor results from 0. 

Before continuing, we briefly document recent attempts to address 
the problem of remote estimation via sensor scheduling. We first ex¬ 
amine offline schemes where sensors are scheduled based on system 
parameters prior to use. Yang et. al. 0 determined that given fixed 
communication constraints, an optimal deterministic offline schedule 
should allocate sensor transmission times as uniformly as possible 
over a finite time horizon. Moreover, Shi et. al. 0 specifically con¬ 
sidered the 2-sensor problem with bandwidth constraints and found 
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that a periodic sensor schedule minimized average error covariance. 
In addition to offline designs, previous work has considered event- 
based designs, where sensor transmissions are scheduled in real time 
based on an occurrence related to a sensor measurement or current 
system parameters. Astrom and Bemhardsson 0 show that for 
certain systems, event based sampling offers better performance than 
periodic sampling. Additionally, Imer et. al. consider a single 
sensor sequential estimation problem where the state is represented 
by an independent identically distributed (i.i.d) process. The authors 
assume communication is limited over a finite horizon and propose a 
stochastic solution. Furthermore, Xu et. al. consider scheduling 
a single, smart sensor which computes and sends a local estimate of 
the state. The authors propose a stochastic event trigger, where the 
rate of transmission is a quadratic function of the difference between 
the state estimate computed at the sensor and the estimate computed 
at the remote estimator. 

While not utilized in 0, IT^ . and (Till , event-based approaches 
can allow the estimator to extract information about the state from 
the absence of a measurement, and thus improve its estimate. For 
instance, Ribeiro et. al. m require the transmission of a single bit 
per observation based on the sign of the innovation and derive an 
approximate minimum mean squared error (MMSE) estimator. Also, 
the authors in fl^ design a threshold scheme on the normalized 
innovation vector to trigger communication to the remote estimator, 
and derive an approximate MMSE estimate. Deterministic schemes 
as discussed by in, in destroy the Gaussian property of the 
innovation process in traditional Kalman filtering, thus rendering the 
closed-form derivation of the exact MMSE estimator computationally 
intractable. Symmetric triggers such as those proposed in Cl and 
0 allow the remote estimator to compute an MMSE estimate. 
Here, the triggers are designed so that a priori and a posteriori 
estimates are identical if a measurement is dropped which implicitly 
requires that the sensor has access to the same information as the 
estimator. However, this is not feasible in the multi-sensor case 
without substantially increasing communication in the network. 

Han et. al. in 0 incorporate a stochastic decision rule, which 
not only allows the remote estimator to use information contained 
in the absense of a measurement, but also maintains the Gaussian 
distribution of the current state. A key advantage of the proposed 
method over most deterministic triggers is that in addition to ob¬ 
taining an exact MMSE estimator, by preserving Gaussianity, 0 
maintains an exact distribution of the state Xk and the estimation 
error et for all time k. Thus, the proposed stochastic event-based 
trigger is useful in scenarios where real time error analysis is critical. 
In this paper, we extend the same stochastic decision rule to the 
multi-sensor case where there exists a unique decision variable for 
each of m sensors. The main contribution of this paper relative to 
0, which considers a binary transmit or drop policy for a single 
trigger, is the derivation of a two-step estimation filter to account for 
multiple independent triggers, a modified optimization problem to 
design each trigger, and a realistic simulation example on data center 
energy management. For this scenario, we also obtain expressions 
for sensor communication rates and upper and lower bounds on the 
error covariance. A preliminary study for this paper was previously 
presented d. Here a three-step recursive filter is proposed which 
computes a state distribution conditioned on all previous information, 
newly received measurements, and the identity of sensors which do 
not transmit sequentially. In this article, we obtain an equivalent two- 
step recursive filter which combines the last two stages, allowing 
us to directly obtain an a posteriori state distribution without any 
intermediary steps. We also extend in by accounting for vector 
sensor measurements with correlated sensor noise as well as through 
our optimization problem and simulation example. 
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The remainder of the paper is organized as follows. Section HU 
formulates the multi-sensor state estimation problem and proposes a 
stochastic event-hased sensor scheduling scheme. Section m intro¬ 
duces a recursive filtering algorithm to obtain the MMSE estimator 
of the state and its error covariance. Section HVl derives results about 
communication rate and estimation performance. Section Ivl proposes 
a semi-definite program to intelligently select trigger parameters. Sec¬ 
tion |Vl] consists of a simulation. A conclusion at the end summarizes 
future work. 

Notation: X' denotes the transpose of matrix X. §+ and §++ 
are the sets of n x n positive semi-definite and positive definite 
matrices. When X € §+, we simply write X > 0 (or X > 0 if 
X G §++)• denotes a Gaussian distribution with mean 

and covariance matrix E. E[-] denotes the expectation, Pr(-) denotes 
the probability of a random event, p(-) denotes the spectral radius of 
a matrix. diag(Xi, • • • , Xs) is the block diagonal matrix with square 
submatrices Xi, ■ • ■ , Xs. 1 and 0 denote vectors with entries 1 and 
0 respectively and In is the identity matrix of size n x n. Finally, 
{A}o is the matrix obtained by deleting all 0 rows from the matrix 
A. 


II. Problem Setup 


We define the following linear system: 

Xk+i = Axk+Wk, t/i*'= i = 1, ■ ■ • ,m. (1) 

Here Xk € M" is the state vector, while is the ith 

of m vector sensor measurements. In addition, Wk € R"' and 


A r (1)/ 

vk = K, 




are mutually uncorrelated Gaussian 


noises with covariances Q > 0 and R > 0, respectively and s = 
To simplify notation, we define yk = , ■ ■ ■ , 

The initial state xo is zero-mean Gaussian random variable with 
covariance matrix Eq > 0, and is uncorrelated with Wk and 
for all k > 0. We assume that {A, C) is detectable where we define 

To reduce the rate of sensor to estimator communication, we 
intelligently transmit a fraction of our sensor measurements. Note 
that we choose to transfer sensor measurements as opposed to local 
estimates. This reduces computation by the sensor as well as possibly 
the size of packets for n > Si. We specify 7 ^“^ G {0, 1} as the 
binary decision variable for sensor i at time k. When 7 ^*^ = 1, 
a transmission occurs while when 7 ^*^ = 0 , no measurement is 
sent. Collecting our decision variables over m sensors, we have 
~ ! ■ ■ ■ 1 Also, suppose at each time k, Ik sensors drop 

their measurements and m — lk sensors transmit their measurements. 
The sensors which transmit have indices pi, • ■ • Define the 

vector of received measurements yj G at time fc by y^ = 

To obtain a MMSE estimator given all previous and current 
measurements, we perform a two-step process. The first step is a 
time update where we obtain the MMSE estimator of Xk given 
the information set up to time k — 1. This is denoted hy Tk-i = 
{ 70 , • • • , 'yk-i,yo, ■ ■■ , yk-i} where I-i = 0. In the second step, 
we update our estimate of Xk, using our previous information set, 
the received measurements at time k, (yk), and the knowledge that 
certain sensors did not transmit a measurement at time k, ( 7 *,). Thus, 
we update using Ik- 

Given the information set, we define the following estimation 
parameters: 


= 'E[xk\Tk-i\, 
Xk = 'R[xk\Tk], 


Pk ='^[{xk-Xk){xk-Xk)\Tk-i\, 

Pk ^'Rlixk - Xk){xk - Xk)'\Ik\- (2) 


Here is an a priori MMSE estimate and Xk is an a posteriori 
MMSE estimate. When all measurements are sent to the estimator, 
computation of Xk and Pk, the error covariance, reduces to the 
standard Kalman filter, where the Gaussian distribution of the state 
allows for a simple recursive filter. As done by 0], to maintain 
the Gaussian distribution of Xk, we consider a stochastic trigger. A 
stochastic trigger takes a measurement y^*^ and computes a function 
^(i) . |Q^ jQ determine the probability sensor i does not 

transmit. While deterministic triggers assign probabilities equal to 1 
or 0 for each measurement, the chosen trigger assigns probabilities in 
[0,1]. To do this, at time k, each sensor i generates an i.i.d. uniform 


random variable over [ 0 , 1 ] and computes 7 ^' 


(») 


(i) 

Tk 


< dHdd) 




(a?) = exp (--otY 




(3) 

Here G S^‘_|_ are trigger parameters and we define Y G S++ 
as y = diag • • • , Note that P{dd = 0|y^“^) has the 

shape of a scaled Gaussian distribution. In the next section, we will 
show this allows the state to remain Gaussian. For the chosen trigger 

we consider stable systems, i.e. p{A) < 1. If the system is unstable, 

(i) 

any sensor i which measures an unstable state will have y^ grow 
unbounded. In this case, by sensor i will always transmit. □ 


HI. MMSE Estimator Design 
I n this section, based on the design of we obtain a closed- 
form solution to the MMSE estimation problem, given recursively by 
the following theorem: 

Theorem 1. Consider remote state estimation with event- 
based scheduler © and define the matrix '^k € R®^® = 

diag( 7 i/sj ■ ■ ■ ■ymism) ft) store the m decision variables. Assume 
f{xo\P-i) ~ A/'(0, Eo) so Xq = 0, Pq = Eq. Then, f{xk\Ik) ~ 
M{xk,Pk) and f{xk\Ik-i) ~ N{x'^,P^) where Xk,x'^ and 
Pk, Pk satisfy the following recursive equations: 

Time update: 

Xk = Axk-i, P^ = APk-iA'Q, (4) 

Measurement update: 


Xk = Xk 

+ PkC'iCPkC + R+{I- ffk)Y-Y\'i’kyk - Cxi), 

(5) 

Pk = Pk - PkC'iCPkC + R+{I- <i)k)Y-^)-^CPk, ( 6 ) 


Proof. To simplify the proof of the theorem, we define the following 
notation which will allow us to distinguish among parameters associ¬ 
ated with sent measurements versus dropped measurements. Suppose 
at time k, there exists Ik sensors ji, - ■ ■ , that do not trigger a 
transmission and m — lk sensors, pi,--- ,Pm-i^, which trigger a 
transmission. We define the matrix Tj, G r(^i=i ''*’0^'’^and Tj, G 
select sensors which transmit and Ak G r(^i=i “^0 
and Ak select sensors which do not transmit as 




Afc = {/ — 


fr ) A J 1 ^ 

u,v otherwise ’ 

{Ak) a|i ^ = . 

10 otherwise 


(7) 


^In (Hi, a closed loop design is considered where a = ]. 

This design can handle unstable systems, but requires estimator to sensor 
communication at each step, which increases communication. As a result, we 
do not consider this approach. 
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We prove Theorem [T] using induction on the distribution 

f{xk\Ik-i) 

Case n = 0: For n = 0, we have Ik-i = 0- Thus, /(a:o|T-i) = 
f{xo) ~ A/'(0, Eo) and the initial conditions holds. 

Case assume for n — ki We assume that f{xk\1k-i) ~ 

m£k,Pk)- 

Case prove for n = k + 1: We first verify the measurement update 
step. 

Measurement Update Step: Consider the joint conditional pdf of Xk 
and KkVk given Xk 

f{xk,f^kyk\Xk) = f{xk,^kyk\Vk,'^k'lk = 1, = 0,Xfe_i) 

~ f {Xk^ ^kyk\Xkyk^ ^k^k ~ —l), 

_ Pr(Afc7;; = 0\xk,yk,Xk-i)f{xk,^^kyk\Xkyk,Xk-i) 
Pr(Afc7fe = 0\Vkyk,Xk-i) 

( 8 ) 

The second equality follows since the knowledge of the values of sent 
measurements Tkyk implies that the decision variables Tk'yk = 1- 
The last equality is derived from Bayes rule. 

By our induction assumption, f{xk, A.kyk,Xkyk) is jointly Gaus¬ 
sian distributed given Xk-i- As a result, the conditional distribu¬ 
tion f{xk, Akyk\Xkyk, Xk-i) is also Gaussian. We first observe 
f{xk,Akyk,Xkyk\Xk-i) has mean [xjl', (AfeCij:)', (PfeCxj:)']' 
and covariance 

Pk PkC'A'k PkC"r'k 

AkCP^ AkiCP^C+ R)A'k AkiCP^C'+ R)r'k . 

PfcCPfe- rk{CP^c' + R)K VkiCP-C+ R)r'^: 

(9) 

Given a joint Gaussian distribution f{xk, Akyk,Xkyk\Xk-i), it is 
easy to compute 

f{xk,Akyk\Xkyk,Xk-i) which is also Gaussian (3]. The condi¬ 
tional means are 

y. = Xk+ Pk{rkCy{rk{CPkC’ + R)r'k)-^rk{yk - Cxi), 

( 10 ) 


Py — AkCXj^ 

+ AkiCPkC + R)r'k{rk{CP^c' + R)r'k)-^rk{yk - Cx^). 

( 11 ) 

Furthermore, the covariance of Xk and Akyk given Tkyk and Xfc-i 


is (fit = 


^xx ^xy 

T' y 

^Xy r-tyy 


where 


S.. = P,- - Pk{TkC)'{Tk{CP^C' + R)T'k)-\TkC)Py 


( 12 ) 


T:yy^Ak{CPkC' + R)A'k- 

AkiCPkC + R)T'k{Tk{CPkC' + R)T'k)-^Tk{CPkC' + PjA^, 

(13) 

E,, = P,-(AfeC)'- 

Pk{TkC)'{Tk{CPkC' + R)T'k)-^Tk{CPkC' + PjA^. (14) 

Now that we have obtained f{xk, Akyk\Tkyk,Xk-i), we also 
observe that 

Pr(Afc 7 fc = 0\xk,yk,Xk-i) = Pr(Afc 7 i; = 0\Akyk) 

= exp (^-^y'kA'kAkYA'kAkyt^ ■ 

(15) 

Using we can thus obtain the joint probability density function 
for the state Xk and the dropped measurements Akyk - That is we have 


f{xk,Akyk\Xk) = Pk ^ exp(-i6»j;), where /3fc € R and 61*, G R are 
defined respectively as 

pk = Pr(A*7fe = 0|r*y*,Xfc_i)\/det(-l>fc)(27r)'‘+sii “H , (16) 



/ 

^xx 


-1 

Xk fJ'X 

-^kVk l^y 


^yx 

'^yy 


^kVk My 


+ {AkVkY AkY Ak{Akyk) ■ (17) 

We now introduce the following Lemma with proof found in the 
appendix. 


Lemma 1. The scalar 0* G R is given by 


6k = 


Xk - Xk 
Akyk — yk 


e: 


Xk - Xk 
Akyk — yk 


=1 , c* G R and 6* G S: 


where Xk G R",i/fe G ^ 
given by 

Xk = x~ + 

Ti—rr! /rr Ti—rr! . 


yk = [^ + AkYA'k] My? ttfe - PyiYyy -f AkY Ak) py. 

( 20 ) 


+ Cfc, 

(18) 


+n 


are 

1 

<H 

1 

), (19) 


0fe = 


Q'xy,k Qyy, k 


( 21 ) 


where 


exx,k = Pk - PkC'icPkC' + R + {i- 'fk)Y-Y^cpp, 

Tilxy,k ~ Yxyi^I ~\~ AkYAkXjyy'j , 

Qyy,k = + A*UA*.] . (22) 


Thus, the joint pdf of our state and unknown measurements are 
given as follows 


f{xk,AkVk\Xk) = exp (-y) 


( 1 

Xk - Xk 

V 2 

^kVk Vk 


Xk - Xk 
Akyk — Vk 


■ (23) 


Since f{xk,Akyk\Xk) is a pdf, its integral normalizes to one which 
implies that f{xk, Akyk\Xk) are jointly Gaussian. Moreover, this 
implies that Xk is conditionally Gaussian given I* with mean Xk 
and covariance P*. Therefore, H) and llUl hold for the measurement 
update step. 

Time Update Step: We have proved /(a:fc|I*) ~ Af{xk, Pk)- By 
the conditional independence of Xk and Wk, we can verify the time 
update step 


f(xk+i\Xk) = f{Axk + Wk\Xk) ~ Af{Axk,APkA' -f Q). (24) 


Thus, 1 ( 4 ) holds. By induction, f{xk\Xk-i) ~ J\f{Xf.,Pf. ). More¬ 
over, from this result, and the proof of the measurement update step, 
f{xk\Xk) ~ M{xk, Pk), which concludes the proof. □ 


Remark 1. The estimation filter can be formulated as a Kalman 
filter with time-varying sensor noise R + {I — ^k)Y~ and inno¬ 
vation ^kyk — Cxf- similarity between the stochastic schedule 
and Kalman filtering allows for computational simplicity and easy 
implementation- 

Remark 2. With an imperfect channel, the estimator will have to 
differentiate between intended packet drops by the sensor due to the 
stochastic trigger and unintended drops due to the channel- If packet 
drops are IID Bernoulli, the state will be distributed according to a 
Gaussian mixture model corresponding to each possible trajectory of 
yk- The resulting distribution however is intractable as k ^ 00 . 
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IV. Performance Analysis 


In proposing an event-based trigger, our goal is to address the trade¬ 
off between estimation performance and power consumed through 
communication by sensor nodes. 

The communication rate £ [0,1] for sensor i can be defined 
as 

A^®' = limsup —^ ^ Ebi*']- (25) 

Knowledge of the communication rate A^®^ of each sensor will 
allow designers to determine the required system bandwidth and 
to estimate the lifetime of each sensor. To obtain an expression 
for the communication rate A*-®^ for each sensor, we first define 

E €§!;:+, £sy+ by 

E = lim Cov(xfc) = AEA' -f Q, 

k —>-oo 

n^®> = lim Covivi"'^) = 


where 7?^®^ = With these results, we now can arrive at 

an expression for the communicate rate of each sensor with proof in 


Theorem 2. Consider a stable linear system o with a stochastic 
event-based sensor schedule given by 13. The communication rate 
A*-®^ for each sensor i = 1, ■■ ■ ,m is given by 


A<^) = 1 - ^ 

^det (7-bn(®)y(®)) 


(26) 


We next verify that the properties established for the expected 
communication rate over several runs, apply to a single sample path, 
the proof of which is found in @]. 


Theorem 3. The following equality almost surely holds. 


lim 

N—^OO 


1 

tTT 


k=0 


(27) 


Furthermore, for any finite integer I > 0, define the event of I 
sequential packed drops over all m sensors Ek,i and the event of I 
sequential packet arrivals over all m sensors Sj, j as follows 

Ek,i = {yk = 0, • • • , 7fe+i-i = 0}, 

^k,i - {7fe = 1, • ■ • , yk-vi-i = !}■ 


Then almost surely j and Ek,i happen infinitely often. 

We next examine the estimation performance by analyzing the 
statistical properties of . 

Theorem 4. Consider a stable system Q with scheduler given by 
id. Let 


gw{X) = AXA' -b Q - AXC'(CXC' + W)~^CXA'. 


1) There exists an M C such that for all k, P^T is uniformly 
bounded above by M. 

2) For any e > 0, there exists an N such that for all k > N, the 
following inequalities hold 

X- el < Pif <X + el, (28) 

where 2L tttid X are the unique solutions X = gn^X) and 
X = gR+Y-''^{X) respectively. 

3) For any e > 0, almost surely for infinitely many k's, we have 
Pjf > X — el and almost surely for infinitely many k's, we 
have Pf < X -\- el. 


The first statement shows that regardless of the choice of 
(communication rate), the error covariance is bounded. The second 
statement obtains upper and lower bounds while the third statement 
shows that during a sample path, P^T will approach these bounds 
infinitely many times, a consequence of Theorem^ where we expect 
long strings of transmissions and drops. 


V. Optimization of Trigger Parameters 
Before we continue, we introduce the following Corollary with 
proof found in the appendix. 

Corollary 1. Define P = X - XC'{CXC' + R + Y-^)-^CX. 

1) For any e > 0, 3 an N such that for all k > N, Pk < P + eL 

2) For any e > 0, almost surely for infinitely many k's, we have 
Pk>P -el 

Thus, it is a worthy goal to design Y^'^ to limit P. We address 
the estimation and communication tradeoff by minimizing the system 
communication rate subject to this bound. 

m 

Problem 1: vj®^ = arg min A^®\ 

subject to P < A. (29) 

0 Here, the matrix A serves as an upper bound on our worse case error 
covariance, thus providing a robust bound on our estimation quality. 
Unfortunately, this formulation deals with a nonconvex minimization 
problem which cannot easily be solved. However, we observe the 
following result. 

Lemma 2. Define f{x) = 1 —(l-|-a ;)“2 and g{x) = 1 —exp(a;)“ 2 . 
Given from ( 125b . H^®^ > 0 and y ^®^ > 0, the following inequality 
holds 

f ^^tr (n(®)y«) j < A°'®‘ < mg 

(30) 



where A°^* is the global minimum of Problem 1. 

Proof. Let u = 'YllLi = tr (n*'®^yi“^). We observe that 

m m 

fiu) < ^ /(Mi) < A°^* < ^ g{ui) <mg(^^ (31) 

2=1 2=1 

The first equality holds for m = 0. The inequality holds since partial 
derivatives of with respect to Ui are greater than or eqml 

to those of /(m). The second and third inequalities are proved in |^. 
Applying lensen’s inequality to g which is concave, we get the last 
inequality. □ 


Since the optimum value of our objective function can be bounded 
by two increasing functions of 'YllLi tr (n*'®^y*'®^), we propose the 
following convex relaxation to Problem 1. 

Problem 2: yj®^ = arg min 
y(*)>o, 

subject to P < A. (32) 

There exist challenges with the constraint since P is only defined 
through X which itself is defined through an implicit function 
gw . The following theorem allows us to obtain an equivalent set 
of constraints and thus formulate the problem as a semi-definite 
program. 


^tr (n«y«), 


2y(t) > 0 is chosen to ensure the problem is feasible for solvers. To ensure 
Y E consider > el where e > 0. 
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Theorem 5. The optimal satisfying Problem 2 can be found by 
solving the following problem. 

m 

Solve: Yj"'> = arg min ^ tr V 

y(‘) >0, 1 = 1 ,,m ' ' 


' Q-^ - S + C'R-^C 

A'Q-^ 

R-^C 


Q-^A 

A'Q-^A + S 
0 


C'R-^ 

0 

Y + R-^ 


> 0 , 


> 0, 5>A"\ 


The proof is found in the appendix. 

VI. Numerical Analysis 



To assess performance, we consider a thermal model for data 
centers, introduced in fl^ . The size of data centers has been growing 
both in number and capacity, resulting in rising energy costs. To 
conserve energy, considers the following thermal model for 
energy control. 


~rpOUt' 

S 


'ksi'^ss - 1 ) ks'iisc ks'^so 


~rpOUt~ 
^ S 

rjiOUt 

C 

= 

ka'^as kai^ac “ 1 ) 


rpOUt 
^ C 

rjiOUt 

. o 


ko'^os ka-^oc kai'^ao “ 1 ). 


rpOUt 


+ Bu, 
(33) 


“rriin" 
^ s 


'Fsa 'Fsc 'Fso’ 


~rjiOUt~ 

s 

rjlin 
^ c 

= 

^cc 4'co 


rjiOUt 

C 

r-pin 

. o 


^oo_ 


rjiOUt 


Here the state a; is a collection of output temperatures of devices 
while the measured values y are the input temperatures of devices 
which require multiple sensors. The subscripts represent different 
nodes under consideration, where ‘s’ corresponds to servers, ‘c’ 
corresponds to air conditioners, and ‘o’ corresponds to other devices. 
The inputs include a reference temperature for the air conditioners, 
power consumed, and temperature of heat sources. gives weight to 
how the temperature output of each node affects the temperature into 
each node and fc is a set of thermal constants. Addressing the trade-off 
between estimation and communication in this example will reduce 
energy expenditures and data storage necessary for thermal control. 

To obtain a model consistent with Q, we linearize the system 
around its stable equilibrium, and assume the inputs remain at or near 
their equilibrium values for all time, a valid assumption during the 
night or backup periods. Furthermore, we sample the system at a rate 
of -^Hz. We consider a system with 16 servers, 3 air conditioners, 
and 1 other device. The matrices Q and R are generated as a product 
of a random matrix with entries uniform from 0 to 1 multiplied by 
its transpose. The matrices are scaled so that the average magnitude 
of error in Wk is 0.1 Kelvin and in Vk is 0.5 Kelvin. In Fig 1, we 
plot the mean squared error in the state estimate as a function of 
the average communication rate, where each data point is obtained 
over a run of 10,000 trials. We consider 3 main designs. We first 
consider a random design where for each sensor at each time step, 
the probability of transmission is Xavg- We also consider a stochastic 
design where each sensor communicates at the same rate, and an 
optimized design from Problem 2. Also shown are upper and lower 
bounds for the un-optimized approach. In Fig 2, we plot the percent 
improvement of the stochastic designs relative to the random design 
in terms of the mean squared error plotted in Fig 1. An un-optimized 
design provides as much as 15% improvement, while the optimized 
design offers as much as 30% improvement. 


Fig. 1. trace(E[Pj. ]) under three scheduling strategies: random, stochastic, 
and stochastic optimized vs Xavg, the communication rate 



Fig. 2. Percent improvement of designed stochastic triggers compared to 
random triggers 


VII. Conclusion 

In this paper we considered a stochastic event trigger for the 
sensor scheduling problem in multi-sensor networked systems. The 
stochastic trigger has inherent advantages over offline triggers which 
can not improve estimates using information contained by the absence 
of a measurement. Moreover, it maintains the Gaussian properties of 
the state, an advantage over previous event triggered approaches. We 
thus could derive a recursive filter to obtain the MMSE estimator 
and error covariance. Additionally, we obtained an expression for 
sensor communication rate as well as asymptotic bounds for our 
error covariance. Finally, we introduced an optimization problem 
that will allow designers to reduce the overall communication rate 
in the system subject to some upper bound on the worst case 
error covariance. Future work consists of considering the stochastic 
trigger in a system with control inputs and incorporating inter-sensor 
cooperation. 
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VHI. Appendix 

Proof. ("Lemma 1) 

We begin by observing the following identities associated with Kk 

and Lfe. 


1) TfeLfc = I , 

2) NkAk = k J sp . ’ 

3) r'j,rfc = im — A'j,Afc = vkfc, 

4) AfcLl, = 0. 

We note that sum of following two quadratic forms can be expressed 
as 


{x — pi)'T,i{x — p,i) + x'T, 2 X 
= a;'(Ei -I- S 2 )a; — 2a;'(Si^i) -|- piSi/n, 

= {x- (El + E 2 )-^Ei/ri)' (El -f E 2 ) {x - (Ei -f E 2 )-^Eipi) 
-t- ^1 (El — Ei(Ei -t- E 2 ) Ei)/ii. (35) 


Setting Y = 


0 

0 


0 

AkYA'k 


, and noting that 


y'kA'kAkY A'kAkVk 


Xk 

/ - 

AkVk 



0 

AkYA'k 


we can directly apply 02} to 03 to obtain 


Xk 

AkVk 


Ok 


Xk - Xk 
AkPk — Vk 



Xk - Xk 
AkPk — Vk 


Ck, 


(36) 


where we have 



Xk 


f^X 

V ^ / 

Vk 


My 


(37) 


Ck 


l-^X 



My 

(38) 


We now attempt to verify ((ID. From Lemma 1 of @], we can directly 
obtain 


— YxyCYyy + AkY ^Tl'^y 

. [Y^yil + AkYA'kYyy)-^]' 


Yxy{I AkYAkYyy) 

[Y-y^ + AkYA'k]-^ ■ 

(39) 


Defining Zfc ^ CPkC + R and Wk = r'kCTkZkW-^Fk and 
substituting equations ([l2ll,([l3l(,([H, we obtain 

Yxx — Yxy (Yyy -|- (AfcYAj,) ) Yx:y 

= P,- - Pk-C'(Vk + Uk)CPk , (40) 

where 

Vk = Wk + {A'k - WkZkA'k) 

X (Afc(Zfc — ZkWkZk)A'k) ^ (Afe — AkZkWk), (41) 
Uk = -{A'k - WkZkA'k) {AkiZk - ZkWkZk)A'k)~" 

X (AfcYA'fc -b {Ak{Zk - ZkWkZk)A'k)~"y ^, 

X {Ak{Zk — ZkWkZk)A'k) ^ {Ak — AkZkWk). (42) 

Since AkZk{Ak — WkZkAk) {Ak{Zk — ZkWkZk)Ak) = I and 
Ak has a unique right inverse equal to A'k we have, 

(A; - WkZkA'k) {Ak{Zk - ZkWkZk)A'ky^ = Z^A'k. (43) 
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Furthermore, we observe that 

{KkZ-;:^h'k){Kk{Zk - z^WkZ,,)Mu) 

= ^kZ^. — AkZf. ^A'l^AkZkWkZkA'^., 

= AkZ^\l - Firk)Zk)K 

— AkZf, '■(/ — F'}J'k)ZkFk{FkZkF'k) ^FkZkAk, 

= AkA)^ — AkZf. Fi^FkZkA)^AkZ)^ Fi^FkZkA/^ 

- AkFiiFkZkF'^y^FtZkAl 
= I. (44) 


Thus, from HD , we obtain 

Vk = Wk + Z^A'^Ak - AkZkWk), 

= Z^A'kAk + (7 - zyA'kAkZk)Wk, 

= zyy^Ak + (7 - zy{i - r'yyzyWk, 

= Z,-iA',Afc + ZyF'yA{F'yF^,Zyy-y^,), 

= Z^A'^Ak + Z^F'^k, 

= zy. (45) 

Moreover, from J43t and l l44b . we have 

Uk = -ZyA'yAkYA'k + AkZ^A'ur^AkZy. (46) 

From the matrix inversion lemma 


Uk + Vk = iCP^C' + R+{I- 
As such, (O is verified. Next from 07} 


Xk 

= (l + ^kY) 


fjk 

\ / 

. My 


7 'ExyiAkY ^A'f. + Yyy) 
[0 {I + YyyAkYA'y-^ 

Therefore, it can be seen that 


My 


(47) 


7 O' 

1 

^xx ^xy 


'0 0 


f-^x 

0 7 

H" 

Y' T 

/-‘xy ^yy_ 


0 AkYA'k 

) 

My 


(48) 


Xk=Xf. + Pf. C'Sk, 


(49) 


where 


Sk = Wk(yk - CxZ) 

- {A'k - WkZkA'y {Ak{Zk - ZkWkZk)K + AkY-’^A'k) 
X {{Ak — AkZkWk)Cxf. + AkZkWkVk)- 

Noting that Wk = WkFkF'k and Aj;rJ. = 0, we have 

Sk = (Wk - {A'k - WkZkA'k){Ak{Zk - ZkWkZk)A'k 
+ AkY-yy-^Ak - AkZkWk)){F'kFkyk - CxZ). 


Utilizing the matrix inversion lemma 


Sk = {Uk + Vk){^kyk-CxZ). 


Thus, 09} and l l20b associated with Xk and yk immediately follow. It 
now remains to verify the expression for Ck in I l20b . We first observe 
by the matrix inversion lemma that 




AkYA'k [0 7] 


= ^k-^k 


(^[0 7]-I>fe ° +Afcy-iA'fc^[0 7]-!-;,. 


Applying O^ . 


Ck = 


f^x 

/ - 

My 



— (E 


0 0 
0 {Yyy+AkY-y'kY 

yy + AkY ^Afc) fj,y. 



f^x 


My 


□ 


Proof. (Corollary 7) To begin we define function h(X, T') : Sj(.+ ^ 
{0,1}“ ^ Sl+ as 

h(X, 4 X - XC' {CXC' + 77 + (7 - diag(T'))y^)”^ CX. 

(50) 

Using the matrix inversion lemma 

h(X, = (^X-i + C' (77 + (7 - diag(T'))y-i)c) . 

This implies h is monotonically increasing in X, and maximized for 
= Os. From Theorem 1, we observe that 

Pk = h{Pk,[YkhY---Yr^i'sJ)- (51) 

By Theorem 4.2, we have that Pff < X + for fc > 7V(e). By the 
monotonicity of h, we obtain 


Pk < h(X + II, 0) 

= (X+77)-^ ^ 

(X + U)C'{C(X + 77)77' + 77 + y^)"^C'(X + 77). 


Moreover, by the continuity of ft in X, for any e > 0 there exists 
7 > 0 such that 


Pk < h(X + U,0) 

< H -XC'iCXC' + R + Y~Y~^C'X + el 
= P + £7, (52) 

for k > N{i) = N{e). We must now show that Pk approaches 
this upper bound infinitely many times. To do this, define function 
h-.Sl+^ Sl+ as 

h{X) ^ (AXA' + Q)- 

(AXA' + Q)C'(C(AXA' + Q)C' + 77 + y ^)-^C(AXA' + Q) 
= h(AXA' + Q,0). (53) 

Note that h is monotonically increasing in X since h is monotonically 
increasing in its first argument and AX A! + Q is monotonically 
increasing in X. Utilizing Proposition 1 of @], we know there exists 
an Z > 0 such that 

h‘{0) >X-XC'{CXC' + R + Y~^)~^CX-eI^ P-el. 

If event Ek,i occurs, then we know that 

Pk+i = h\Pk)>h\G)>P-eI. (54) 

By Theorem 3, the event Ek,i almost surely occurs infinitely often 

and thus the result holds. □ 


Proof. (Theorem 5) We first assert that the following two statements 
are equivalent. 

1) P < A, 

2) There exists 0 < U < A such that U < h(U). 

The first statement implies the second by taking U = P. Noting 
the monotonicity of h and the convergence of hA to the fixed point 
P, the second statement implies the first by repeatedly applying h. 
Take S = U~^. Then by the matrix inversion lemma, the following 
statements are equivalent. 

1) P < A, 































2) There exists 5" > A ^ such that 0 < (AS" ^ A' + Q) ^ + 

c'iR + Y-'^y^c - S. 

By the matrix inversion lemma, 

{AS-^A' + Q)-i = g-i - Q-^A{S + A'Q-^A)-^A'Q-\ 


Thus, using Schur’s condition for positive definiteness we have that 
the following statements are equivalent. 

1) P < A, 

2) There exists S > A“^ such that 

'Q-^+C'{R + Y-^)-^C - S Q-^A 1 

A'Q-i A'Q-'^A + S\- ' 

(55) 


A'g-^A + 5 > 0. 


Since S' > 0, the latter statement immediately holds. Now by the 
matrix inversion lemma, 155) is equivalent to 

'Q-^ +C'R-^C - S Q-^A 

A'Q-i A'Q-^A + S 

{R~^+ Y)-^ [R-^C 0 ]> 0 . 


Thus, by Schur’s condition for positive definiteness we have that the 
following statements are equivalent 

1) P < A, 

2) There exists S > A~^ such that 


Q-^ + C'R-^C - S 

Q-^A 

C'R-^ ' 


A'Q-i 

A'Q-^A + S 

0 

>0, 

R-^C 

0 

R-^ +Y 




R~^ + Y > 0. 


Note that the latter statement is given for free since P ^ > 0. The 
theorem follows immediately. □ 










